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ABSTRACT 

This  paper  describes  a  canonical  procedure  to  approximate  an  arbitrary 
family  of  C  vector  fields  {g^,...,gm)  on  S?3  with  vector  fields 
9j»...#gm  on  the  same  space  R*3  Which  generate  a  nilpotent  Lie  algebra. 

Each  g^  can  be  either  obtained  from  a  Taylor  expansion  of  the  input-output 
map  for  the  control  system 

m 

(*)  x  -  l  g. (x)u  ,  x<0)  -  0  , 

i-1  1  1 

or  computed  directly  as  an  asymptotic  limit  of  the  corresponding  vector 

field  g^.  A  useful  consequence  is  that  every  control  system  of  the  form  (*) 

00 

can  locally  be  regarded  as  an  arbitrarily  small  C  perturbation  of  a 
nilpotent  system  on  the  same  state  space,  up  to  a  suitable  linear  rescaling  of 
coordinates . 

AMS  (MOS)  Subject  Classifications:  93b10,  93C10 

Key  Words:  Dilation,  nilpotent  family  of  vector  fields 

Work  Unit  Number  5  (Optimization  and  Large  Scale  Systems) 
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SIGNIFICANCE  AND  EXPLANATION 


If  a  family  of  vector  fields  generates  a  nilpotent  Lie 

I  m 

algebra,  then  the  response  of  the  control  system 


m 

x(t)  -  l  g  <x)u. <t),  x( 0 )  **  0  e  S*5 
i-1  1  1 


can  be  written  out  explicitly  in  terms  of  integrals  of  the  controls  u^.  This 
and  other  nice  consequences  make  nilpotency  a  highly  desirable  property  from 

the  point  of  view  of  mathematical  analysis. 

. 

-  The  present-- paper-  describes  a  canonical  method  to  locally  approximate  any 

►  p  !•  "S  ■  »  ^ 

family  of  vector  fields  by  one  which  generates  a  nilpotent  Lie 

algebra.  This  is  particularly  useful  in  control  theory,  because  it  shows  that 
an  arbitrary  control  system  can  locally  be  obtained  from  a  small  perturbation 
of  a  nilpotent  system,  by  a  suitable  rescaling  of  the  coordinates. 

f 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


LOCAL  ASYMPTOTIC  APPROXIMATION  OF  NONLINEAR  CONTROL  SYSTEMS 


1.  INTRODUCTION 


Alberto  Bressan* 


This  paper  Is  concerned  with  an  autonomous  nonlinear  control  system  of  the  fora 


x(t)  “  l  g,  (x(t) )u  (t),  x(0)  -  0  e  Rd  , 

i-1  1  1 


(1.1) 


where  the  g^  are  C  ,  globally  bounded  vector  fields  on  and 

u(>)  »  (u1 .u^ )  e  L1 ( to," ) iEB) .  The  system  (1.1)  generates  a  smooth  input-output 
map  4  s  u(*)  ♦  x(u,*)  from  L*(  tO,»)rn")  into  C°(  [0,**)*Rd).  Explicitly  computable 
approximations  of  4  are  of  primary  importance  in  the  local  study  of  (1.1).  The  p-th 
order  Taylor  expansion  T®4  of  4  about  the  null  control  was  studied  in  [1.3].  In  [8]  it 
is  shown  how  the  trajactoriea  of  (1.1)  can  be  locally  approximated  by  means  of  an 
additional  (nllpotant)  system,  say 


x(t)  »  l  f  (x(t))u  (t),  x(0)  -  0  e  tfi’  . 

i»1  1  1 


(1.2) 


For  certain  applications,  such  as  the  coeiputation  of  a  local  time-optimal  feedback  [2], 
both  approaches  seem  unsatisfactory.  In  general,  1^*4  is  merely  a  sum  of  multilinear 
integral  mappings  from  L1  into  C°  and  does  not  arise  as  the  exact  input-output  map  of 
any  control  system.  Bis  optimality  of  a  given  control  u  under  Tp4  cannot  therefore  be 
tested  by  the  Maximum  Principle.  On  the  other  hand,  the  vector  fields  f j,  in  (1.2) 
represent  a  lifting  of  the  g^  (i  ”  1,...,m)  in  a  usually  higher  dimensional  space  R*5  . 
It  is  not  always  possible  to  determine  a  local  property  of  (1.1)  by  studying  (1.2),  since 
the  two  systems  live  on  different  spaces. 
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Using  a  singular  perturbation  technique,  in  this  paper  we  derive  a  canonical  procedure 
to  approxiaate  the  vector  fields  g^  with  a  nilpotent  family  {  g  ^ , . . . , g^}  of  vector 
fields  defined  on  the  sane  state  space  fP.  The  main  construction  is  as  follows.  For 
snail  e  >  0,  the  restriction  of  the  map  i  to  controls  with  support  inside  [0,e]  can 
be  rescaled  to  a  nap  fron  L1) [0, 1] ,rf“)  into  C°( [0, 1] itf1) .  This  is  achieved  by 

neans  of  the  transformation  t  ♦  e  t  of  the  time  variable  and  by  letting  a  suitable 
faaily  of  dilations  4^  act  on  the  space  variables.  For  each  e  >  0,  $E  is  thus  the 
input-output  map  generated  by  the  control  system 

B 

x(t)  -  l  g^(x(t))u. (t),  x(0)  -  0  e  Rd  ,  (1.3) 

i-1 

where  the  vector  fields  gE  are  obtained  from  the  g^  after  a  rescaling  of  coordinates. 

As  e  ♦  0,  there  exist  a  nap  ♦  and  vector  fields  g^  on  such  that  ♦  ♦  and 

gE  ♦  g^  together  with  all  derivatives,  uniformly  on  bounded  sets.  ♦  is  then  the  input- 
output  nap  corresponding  to  the  control  system 

m 

x(t)  »  l  g  (x(t))u. (t),  x(0)  »  C  .  (1.4) 

1-1  1 

The  vector  fields  g^  have  polynomial  coefficients  ani!  are  invariant  under  a  1-parameter 
group  of  transformations.  Moreover,  without  any  assumption  on  the  Lie  algebra  generated 
by  gt,.*.,gm,  it  turns  out  that  Lie  (o(,...,g  )  is  always  nilpotent,  so  that  the 
solutions  of  (1.4)  can  be  written  in  closed  form  [7).  This  remarkable  feature  makes  the 
system  (1.4)  an  attractive  object  for  a  detailed  mathematical  analysis.  Indeed,  since 
(1.4)  is  obtained  as  a  uniform  limit  of  (1.3)  as  e  ♦  0,  any  property  of  (1.4)  which  is 
retained  under  small  C*  perturbations  [4]  yields  a  local  property  of  (1.1). 
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2.  PRELIMINARIES 


In  the  following,  we  use  I* I  for  the  euclidean  norm  on  Rd  and  I • I  for  the  norm 

In  Banach  spaces.  The  closed  ball  centered  at  x  with  radius  r  is  denoted  B(x,r), 

while  0  indicates  tensor  product.  Given  two  Banach  spaces  E  and  F,  k  >  0,  we  denote 

by  L  (E;F)  the  space  of  continuous  k-linear  mappings  A  from  0E-E0E0...0E  (k 

k 

times)  into  F  with  the  operator  norm 

IAI  »  sup{l  A(v1 , . ..  >vk)lp!  <  1,  i  -  1,...,k)  . 

If  4>  s  E  ♦  F  is  a  smooth  mapping,  its  k-th  Frfichet  derivative  a  a  point  u  e  E  is 
D%(u)  e  L*(E;F).  The  k-th  order  Taylor  expansion  of  ♦  at  the  origin  is  then 

k  k  -  . 

T*(u>  •  1  5T  D>(0)  •  u1^  , 
y  J-0  3 

where  u^  —  (u,u,...,u)  e  0  E.  The  sane  notation  T*g  is  used  for  the  k-th  order  Taylor 
expansion  of  a  vector  field  g  on  Hr  at  the  origin.  It  was  shown  in  [1]  that  the  n-th 
order  Taylor  expansion  t"$  of  the  input-output  map  ♦  generated  by  (1.1)  about  the  null 
control  can  be  obtained  by  computing  the  n-th  Picard  iterate  for  the  approximate  system 

m 

x(t)  -  ),  Tn-1gi(x(t))u1(t),  x( 0 )  -  0  ,  (2.1) 

i-1 

discarding  the  terms  of  order  >n.  The  multilinear  integral  mappings 
U  :  L1  (<),«•)  *  C°(0,*)  arising  from  the  above  procedure  were  called  in  [1]  integral 
monomials.  If  u  is  k-linear,  we  say  that  y  has  order  k.  The  first  Picard  iterate  for 
(2.1)  is 

m  t 

Pl<u,t)  -  l  j  gi(0)u1(s)ds  • 
i-1  0 

In  general,  if  the  k-th  Picard  iterate  for  (2.1)  is  a  sum  of  integral  monomials,  say 

N(k) 

Pt (u,t)  ■  l  V,  (u,t)  , 

*  i-1  1 

then  P^+1  can  again  be  written  as  a  sum  of  terms  of  the  type 


-3- 


(2.2) 


V(u,t) 


t  .  i 

J  77  D-V  (0)  •  (u  (u,s), 
0  31  1  A1 


,U  <u,s))u  (s)ds 

3 


where  0  <  j  <  n,  1  <  i  <  m.  If  uA  has  ordar  (t  «  1,...,j),  than  v  in  (2.2)  has 
ordar  +  ...  +  vj  ♦  t«  The  abova  construction  canonically  determines  an  Increasing 
sequence  of  subspaees  Zp  C  sf3,  namely  Zg  «  {0},  Z,  »  span{g^(0)i  i  "  1,...»m)  and 
Inductively 


Zp  m  spantD^g^O)  •  (v1>....Vj)>  1  <  i  <  m,  0<j<p, 

Vj  e  zp^.  P1  ♦  ...  +  Pj  <  p)  •  l2’3) 

Comparing  (2.3)  with  (2.2),  it  is  clear  that 

TP*(u)(t)  6  Zp  (2.4) 

for  all  p  >  0,  u  e  I1,  t  >  0.  Indeed,  Zp  is  precisely  the  subspace  of  tP  spanned  by 
the  coefficients  of  the  integral  monomials  of  order  <  p  in  the  Taylor  expansion  of  <|> 
about  the  null  control. 


i 


3.  A  CLASS  OP  NILPOTENT  LIE  ALGEBRAS 


All  of  the  vector  fields  arising  from  the  asymptotic  limit  procedure  considered  In 
this  paper  generate  a  special  type  of  Lie  algebras  which  we  now  briefly  describe.  Por  a 

CD 

different  approach,  based  on  the  dual  action  of  vector  fields  as  derivations  on  C 
functions,  see  16].  Given  an  orthogonal  decomposition  Rd  "  W1  A  ...  A  W-,  let  be 

the  canonical  projection  of  iP  and  the  differentiation  w.r.t.  the  p-th  component 

of  x.  Define  L  (Lr,  for  r  >  0)  to  be  the  set  of  all  smooth  vector  fields 
f  •  (f1(...,f-)  on  iP  such  that 


whenever  k>0,  1  <  q  <  p  and  pj  +  ...  +■  p^  *  q  (pj  +  ...  +  pJc>q-r).  For  *  > 
Lr  contains  only  the  null  vector  field.  From  the  above  definitions  it  follows 

LEMMA  1.  L  is  a  finite  dimensional  nilpotent  Lie  algebra  of  vector  fields  with 
polynomial  components,  with  the  usual  bracket  operation  [f,g]  “  (Dg)  •  f  -  (Df)  •  g, 
and  L  «  Lq  D  L1  D  . . .  2  is  a  decreasing  sequence  of  ideals ,  indeed 


P> 


(Lr,L.]  £  Lr+.+1  (3.1) 

PROOF.  If  f  e  L,  any  p-th  derivative  of  f  vanishes  identically.  Hence  L 
contains  only  vector  fields  whose  components  are  polynomials  of  degree  less  than  p.  In 
particular,  L  is  finite  dimensional.  To  prove  (3.1),  let  f  e  Lr,  g  e  Lg ,  k  >  0, 
q,p1,...,pk  6  {1,...,p}  with  p1  +  ...  +  p^  >  q-r-s-1.  Denoting  fq  -  x^f,  gq  -  x^g  the 
q-th  components  of  f  and  g,  one  has 


9x 


3x 


x  I  (Df ) 

q 


q] 


3x 


if  (r-  ' 

1-1  3xi  q 


)  *  9*] 


(3.2) 


Notice  that  the  right-hand  side  of  (3.2)  can  be  written  as  a  sum  of  terms  having  the  form 


In — 

po(1) 


,h+l 


^o(h) 


f,5  •  W: 


,k-h 


Po(h-M) 


- 

Po(k) 


(3.3) 


-5- 


where 


0  <  h  <  k 


1  <  i  <  p  and  0  ia  a  permutation  of  the  aet  ( 1 


k) .  if 


pc(1)  +  *•*  +  po(h>  +  1  *  q  "  r  ’  tha  firat  term  in  the  tanaor  product 
identically.  Otherwiae  }  *  •••  +  Pff(h)  *  q*r-i-1,  hence 

po(h+1)  +  *  pc(k>  *  -  1  -  i 

and  the  aecond  factor  in  (3.3)  vaniahea.  Therefore  (Df)  •  g  @  Lr+J+  y . 
(Dg)  •  f  e  l<r+8+')#  hence  (3.1)  holde/  proving  the  nilpotency  of  L. 


(3.3)  vaniahea 


Similarly 
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4.  STATEMENT  OF  THE  MAIN  RESULT 


With  the  increasing  sequence  of  subspaces  Z  c  ***  defined  at  (2.3)  one  can  associate 


an  orthogonal  decomposition 


IT  -  W,  *  ...  ®  W- 
1  P 


as  follows.  Fix  any  p  >  1.  For  1  <  p  <  p  let  Wp  be  the  orthogonal  complement  of 

Zp_1  in  Zp ,  and  let  W-  be  the  orthogonal  complement  of  zp-t  in  Thin  clearly 

yields  (4.1).  The  canonical  projection  of  jP  onto  Wp  is  denoted  x^.  In  addition  to 

the  input-output  map  i  :  u(  • )  ♦  x(u,* )  generated  by  (1.1),  for  0  <  e  <  1  we  can  now 

G  G 

define  the  rescaled  maps  ♦  :  ui • )  ♦  x  (u,»)  by  setting 


xE(u,t)  =  \  E  px  (x(eu,t)) 

p-1  p 


By  direct  computation  one  checks  that 


x  (u,t)  »  l  g  (x  (u,t))u.(t),  x  (u,0 )  «  0 
i-1  1 


gj(x)  -  l  e1  px  [g  (  l  (x))]  . 

p«1  v  j-1  1 


Our  major  interest  is  in  the  behavior  of  the  rescaled  system  (4.3)  as  e  ♦  0. 

THEOREM.  Let  g^  ( 1  “  1 , . . . ,m )  be  C  ,  globally  bounded  vector  field#  on  tfl,  and 
let  ♦E,xe,g^  be  defined  by  (4,2),  (4.4),  corresponding  to  the  decoeiposltion  (4.1) 
obtained  from  (2.3).  Then  as  e  ♦  0 

1)  converges  to  ♦  :  u(*)  ♦  x(u,*),  defined  by 


♦<u)(t>  -  \  x  [T*>(u)(t)]  . 


)c  C  )( v 

More  precisely,  for  all  k  >  0,  D  *  tends  to  D  4  uniformly  on  bounded  subset#  of 


l’uo.-wW*). 
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ii)  For  all  1  -  g*  convergea  to  g1#  defined  by 


P  P-1  , 

Ji<x)  “  l  *01  l  ^  Tf  "  n 

1  p-i  p  j-o  oer<P,j)  3 


—  D3g,(0)(icg{1)<x>, 


,*0(j,<x»]  ,  (4.6) 


f(p,j)  being  the  set  of  all  nape  0  i  for  which 

0(1)  ...  +  0(j)  -  p  -  1.  For  all  *  >  o,  DKgj|  tends  to  0  gt  uniformly  on  bounded 

subsets  of  K?3 . 

iii)  por  all  u(*)  e  t1,  the  trajectory  t  ♦  x(u,t)  -  7(u)(t)  is  the  solution  of 


x(t)  •  l  g.  (x(t) )u. (t) »  x(0 )  “  0  . 

i-1  1 

iv)  L.'.e  {g,,...,g  }  is  nilpotent. 

1  w 


('.  7) 
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5.  PROOF  OF  THE  THEOREM 

Fix  any  p  #  and  write 

*  (♦E(u)(t))  «  £  px  ($(cu)(t)) 

P  P 

■  £-px  (Tp*(eu)(t))  +  e"pir  (n(eu)(t)) 

P  P 

where  n  "  *  -  Tp^,  the  remainder  of  the  p-th  order  Taylor  expansion  of  ♦  about  the  null 
control,  is  a  C  map  from  b'  ( (0,» ) ;  R®)  into  C([0,<»)jRd)  with  D^n(O)  «  0  for 
0  <  j  <  p.  By  (2.4),  x^T1*  %  =  0.  Therefore,  11^4  is  a  homogeneous  p-linear  mapping, 
and 

e'Px  1?4(cu)  -  x  tP*(u)  V  £  >  0  .  (5.1) 

P  P 

__  lr  £ 

To  prove  i)  it  now  suffices  to  show  that,  for  all  k  >  0,  e  PD *(x  •  n  )(u)  converges  to 

zero  uniformly  as  u  ranges  on  bounded  subsets  of  l',  nE  being  the  map  u  *  n(cu).  The 

assumptions  on  n  imply  the  existence  of  a  constant  C  >  0  such  that 

IDkn(u)l  <  C  min{lulP_k+\l} 

whenever  lul  <  C  ' .  If  U  is  a  bounded  subset  of  L1,  choose  £Q  >  0  so  small  that 

e0°  -  If  u  e  0  and  0  <  £  £  e0#  then 

i£*pDk(x  •  ne)  (u)l  <  £k-pIDkn(eu)l  <  C£k_p*min{ep'k+1iul ,  1)  .  (5.2) 

P 

In  both  cases  k  >  p  or  k  <  p,  as  £  *  0  the  right-hand  side  of  (5.2)  converges  to  zero 
uniformly  on  U.  This  establishes  i).  The  proof  of  ii)  is  -.imilar:  fix  i  e  {1,...,m} 
and  p  e  (l,...,p).  Using  a  Taylor  expansion  of  g^  of  order  p  -  1  with  remainder 
h(»),  from  (4.6)  one  obtains 


X  g£(x)  -  w  [e’-P  l  -L  l  Djg  (0)(£°(1) 

p  1  p  j-o  jl  oerj  1 


\J(i),x) . EO(j)no(j)(x)>] 


+  x  e1  p  h(  )  £*x  (x) 

p  t-1  * 


(5.3) 


where  is  the  set  of  all  maps  o  from  {1,...,-'}  into  {1,...,p}  and  h  is  a  C 

vector  field  on  R^  with  D^h(0)  -  0  for  j  “  1,...,p  -  1.  The  generic  term  in  the  first 
summation  in  (5.3)  is  a  homogeneous  j-linear  map  from  Rd  into  Wp  of  the  form 
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V(e*x>  ■  Vir  e1'pD391(0)(e<,<1,«<y(1)(x) . ta{lKgii)U))] 


,  i-p-s.  . 

“  *p[-j7  e  DJ,1(0)OraM)<x) . *0(J)<x>>] 


1-p+S 

with  Sg  -  d(1)  +  ...  +  0<j).  If  Sg  >  p,  than  11m  ?<e,x)  -  lim  e  °T(1,x)  -  0  and 

e*0  c*0 

tha  aan a  hold*  for  all  derivatives  of  1  w.r.t.  x,  uniformly  on  bounded  aubaata  of 
F^.  If  S0  ■  p  -  1,  then  F(e,x)  -  F(1,x)  doaa  not  depend  on  e.  If  s  <  p  -  1,  since 
»„,,,<*)  «  z<j(t)  th“  definition  of  implies  T(e,x)  e  Zp_i»  hence 

*  T(e.x)  -  0.  To  prove  ii)  it  now  suffices  to  show  that  the  second  summand  on  the 
right-hand  side  of  (5.3)  converges  to  zero  as  e  ♦  0.  For  every  k  >  0,  the  aseumptions 
on  h(*)  imply  the  existence  of  a  constant  C  >  0  for  which 

|Dkh(x)|  <  C*min{  |x|P  lt,  1}  (5.4) 

whenever  |x|  <  C1.  Let  V  be  a  bounded  subset  of  fP,  _and  choose  eQ  e  (0,1)  so  saill 
that  eQV  _  B(0,C_1).  For  0  <  e  <  eQ,  set  he(x)  -  h(  £  ep*  (x)j.  From  (5.4)  it 

p-1  p 

follows  the  estimate 


e1_p|Dkhe<x)|  -  e,-psup{|Dkhc(x)*ytkJ  |- lyl"k>  y  e  Rd\{o}} 

<  e1  Psupf|Dkh(  l  epit  (x)W  £  tPr  (y)j  [kl  I*  lyl  k»  y  *  o( 

p-1  p  p-1  p 

<  e1  psup{c*min{ |ex|p_k, 1)* ley |k* |y|  k>  y  *  0} 

<  ce1""p+kmin{ep-k|x|p"k,1}  .  (5.5) 

In  both  cases  k  >  p  or  k  <  p,  the  right-hand  side  of  (5.5)  converges  to  zero  uniformly 
on  V,  as  e  +  0.  Since  the  above  holds  for  every  i  and  p.  ii)  is  proved.  Using  i) 
and  ii).  iii)  follows  from  (4.3),  letting  e  +  0  in  the  equality 

,  t  m  ,  . 

x  (u.t)  mJ  l  g  (x  (u.s) )u(s)ds  . 

0  i-1  1 

Finally,  we  check  that  each  g^  belongs  to  the  Lie  algebra  L  defined  in  §3, 
corresponding  to  the  decomposition  (4.1).  Let  k  >  0,  q.p^,...,^  e  {1,...,p}, 


Pi  +  . . .  +  Py  >  q.  By  (4.6),  the  q-th  component  of  g^  is  a  sum  of  terms  of  the  form 


*<x)  -  4-  IT  Ib’g  (0)  (v  (x>,...,«  (X))] 

j*  <!  i  q1 

with  qj  +  ...  t  qj  »  q  •  1.  We  therefore  have 


3k* 


3x 


(x) 


0 


for  all  x  e  sf5  unless  there  exists  an  injection  a  of  the  set  into  the  set 

such  that  fjo j  «  pt  for  all  l  »  1,...,k.  But  the  existence  of  such  an 
injection  would  contradict  the  assumption  p(  ♦  ...  ♦  py  >  q.  Therefore  g^  e  L  for  all 
i  e  By  Lemma  1.  Lie  (g }  is  a  subalgebra  of  a  nllpotent  Lie  algebra. 


This  completes  the  proof. 


6.  CONCLUDING  REMARKS 


The  above  raaulta  also  apply  to  control  ayatana 


a 

x(t)  -  l  X^xttMu^t),  x(0)  -  C  (6.1) 

1-1 

on  a  d-dimensional  manifold  M.  Indeed,  if  tha  Lla  algebra  generated  by  has 


full  rank  at  £ .  then  aaong  iterated  brackets  of  the  X^'s  one  can  choose  d  vector 
fields,  say  Yj,...,Yd,  which  are  linearly  independent  at  K  •  The  nap 

8  i  (sv...,sd)  *  (exp  e,Y,)  •  ...  •  (exp  adYd)(C>  (6.2) 

thus  defines  a  local  chart  of  a  neighborhood  of  €  in  M  (81.  The  asymptotic  expansion 
considered  in  {4  can  now  be  performed  for  the  system  on  R*3  corresponding  to  (6.1)  in  this 
canonical  chart. 

The  present  approximation  technique  enables  one  to  study  -  local  problem  concerning 
(1.1)  by  first  solving  the  corresponding  problem  for  the  system  (4.7)  and  by  then  proving 
that  the  structure  of  the  solutions  is  ’stable”  under  suitably  small  C  perturbations  of 
the  vector  fields  g^ •  This  approach  is  adopted  in  [2]  to  study  the  local  time-optimal 
stabilizing  feedback  for  a  generic  three-dimensional  nonlinear  system  with  scalar  control. 

The  analysis  of  the  control  system  (4.7)  takes  advantage  from  the  explicit 
representation  (4.5)  of  trajectories  in  terms  of  integrals  of  the  controls.  One  can  also 
make  use  of  the  special  properties  of  the  vector  fields  gi  in  connection  with  the  family 
of  dilations  :  Rd  ♦  Rd, 

P  „ 

«,(x>  -  i  ep*  (x> 

*  p-i  p 

(see  (6)  ).  Some  useful  consequences  deserve  mention.  Fix  a  compact  convex  set  fl  C  Rm 
and  define  the  admissible  set  of  controls 

v  -  {u( • )  e  l’( (0,-)>Rm)»  u( t )  e  n  v  t  >  o>  . 

Given  a  control  u  e  U,  for  5  >  0  define  u^(t)  -  u(£t).  Call  R(t)  the  set  reachable 
at  time  t  by  trajectories  t  ♦  x(u,t>  of  (4.7)  with  controls  in  U.  With  these 
conventions  we  have 
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*&}&"+■ 


LEMMA  2.  Let  t  >  0,  t  >  0,  i  e  u  e  U.  Then 


gi<6  (x))  ■  i  !r  t  (g^x))  , 

p-1  P 


Cgi(*?(x))  -  S^fg^x))  , 


x(fu^,t)  -  x(u,£t)  -  S^(x(u^,t))  , 


R<0  -  i?R{1) 


PROOF.  By  (4.6),  ( 5 ( x ) )  la  a  aua  of  terms  having  the  form 

T(VX))  "Vi?  DV°>  •  (\»(i)(Vx,) . 

with  0(1)  +  ...  ♦  o(j)  ”  p  -  1.  Therefore  T(5^(x>)  *  E^Ntx).  Thus  yields  (7.1). 
Multiplying  both  aides  of  (7.1)  by  £ ,  one  gets  (7.2).  The  first  equality  in  (7.3)  holds 
simply  because  (4.7)  is  an  autonomous  system,  linear  in  u.  Moreover 


c(u,Ct)  -  «.(x(u,,t))  -  l  J  t5g.  <*(u,tt) )  -  6_(g  (x(u,  ,t)))]dt 
*  *  i»t  0  1  tit 


Gronwall’s  lemma  and  (7.2)  now  yield  the  second  equality  in  (7.3),  from  which  (7.4)  follows 
by  setting  t  »  1. 

Notice  that  a  control  u  e  O  is  time-optimal  on  t0,£)  iff  uE  is  time-optimal  on 
(0,1).  if  (4.7)  admits  a  regular  time-optimal  feedback  up  and  if  t  ♦  u(t)  is  optimal 
on  [0,T],  then  for  t  e  (0,T) 

up(x(u,Ct))  ■  u(Ct>  *  u^tt)  »  Up(x(u^,t))  ■  up(4^_1(x(u,tt)))  . 

Therefore  we  expect  uv  to  be  invariant  under  dilations. 


-'•ftfettfe**. 


7 .  Ml  EXAMPLE 


Consider  the  two-dimensional  system 
d 

~  (x,y)  ■  (u1  +  sin(x  ♦  y)u2,  (1  -  cos(x  *  y))u2)  .  (7. 1) 

(x(0),y(0>>  -  (0,0)  , 

Computing  the  third  Picard  iterate  for  the  reduced  system 

<x,y)  -  (u,  +  (x  ♦  y)uJt  \  (x  ♦  y)2u2),  (x(0),y(0))  -  (0,0) 
and  discarding  term*  of  order  >  3,  one  obtaina  the  third  order  Taylor  expansion  of  the 
input-output  map  ♦  generated  by  (7.1)t 


3  *  t  a 

T  d(u)(t)  -  ( I  u1(a)d a  ♦  }  ( J  u, (o)do)u,(s)dn 

0  0  0 

t  o,  a. 

*00  4 

,  t  a  _ 

r  J  (J  u  (o)do)  u(a)de)  . 

*  0  0  1  2 

In  thin  case  l,  ■  z2  «  »,  ■  (|*,0)ix  1  «),  Z3  -  R2  , 

W2  -  f(0,0)},  M3  -  { (0,y) iy  e  r}  . 

The  'homogenised'1  expansion  (4.5)  la 

-  t  «  t  a  a 

♦  (u)(t)  »  [J  u1<a)ds,  j  J  (J  u, (o )do j  u2(a)ds) 
0  *  0  0 


which  exactly  represents  the  response  of 

§JT  ( x,y)  -  (u1,1^x2u2),  (x(0),y(0>)  -  (0,0)  .  (7.2) 

If  the  aet  of  admissible  controls  is 

0  -  {u  -  (u^Uj)  e  L1((0,-)iR2)»  I u, ( t )  1  <  1,  u2(t)  -  1,  V  t  >  0}  , 
the  reechable  eeta  for  (7.2)  are 

R(t)  -  { (x,y)  C  R2|  ^  I x | 3  <  y  <  i  ( +  t2  |x|  *  t |x| 2  -  |x|3) 
and  the  relations 
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hold  for  every 


with 


Notice  that  u_ 


<x,y)  e  R(t)  iff  Ux.S'y)  e  R(Ct) 
t,C  >0.  A  time-optimal  feedback  for  (7.2)  is 

Up(x,y)  -  (1,1)  on  A+,  up(x,y)  -  (-1,1)  On  A~  , 

A+  »  {(x,y)j  x  >  0,  y  »  •g  x3}  U  {(x,y)>  x  <  0,  y  >  ~  |x| 
A  »  {(x,y)>  (-x,y)  e  A+}  . 

is  invariant  w.r.t.  all  dilations  5,  :  (x,y)  *  (Cx,C3y) 
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